Icosahedral Ti-Zr-Ni: A groundstate quasicrystal? 
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The first complete ab initio zero-temperature ternary phase diagram is constructed from the 
calculated energies of the elemental, binary and ternary Ti-Zr-Ni phases; for this, the icosahedral 
i-TiZrNi quasicrystal phase is approximated by periodic structures of up to 123 atoms/unit cell, 
based on a decorated-tiling model. The approximant structures containing the 45-atom Bergman 
cluster were nearly degenerate, and stable against the competing binary phases. It is speculated 
that i-TiZrNi may be a ground state quasicrystal, as it is experimentally the low-temperature phase 
for its composition. 

PACS numbers: 61.44.Br, 71.15.Nc, 81.30.Bx 



Thermodynamically stable, long-range ordered qua- 
sicrystals of icosahedral symmetry are known in both 
of the main structural classes of quasicrystal, the Al- 
transition metal class (e.g. i-AlPdMn) and the Frank- 
Kasper class @ (e.g. i-ZnMgY). The Ti-based quasicrys- 
tals such as i-TiCrSiO and the thermally stable i- 
TiZrNi [|[ £| || fall into either respective class. They 
are of technical interest due to their high melting point 
and capacity for hydrogen absorption, but arc less stud- 
ied because the crystallites are small and the coherence 
length is only 35 nm. 

In experiments on the TiZrNi alloy, the crystal struc- 
ture VF-TiZrNi is stable at high temperatures, but upon 
cooling it undergoes a reversible phase transition at 
570° C to i-TiZrNi, showing that i-TiZrNi is lower in en- 
ergy than PF-TiZrNi || |. In fact IF-TiZrNi is a "peri- 
odic approximant" of i-TiZrNi, meaning that the unit cell 
is identical to a fragment of the icosahedral. Long-time 
anneals (up to one month) at 500° C gave no indication 
that the quasicrystal transforms to some other phase. 
The situation here contrasts with the Al-transition metal 
class, in which the analogous crystal a-AlMnSi (known 
as "1/1 approximant"), is lower in energy than the qua- 
sicrystal of identical composition. 

In this work, we address the possibility that the i- 
TiZrNi quasicrystal (or a very large unit cell crystal of 
nearly identical structure) is in fact a ground state. We 
compute the ternary ground state phase diagram, using 
a database of ab initio total energies of 28 elemental, 
binary and ternary crystalline phases in the Ti-Zr-Ni al- 
loy system. The total energy for the quasicrystal phase 
is determined from six different periodic approximants 
of i-TiZrNi. This calculation depends on our previously 
reported structure model which was formulated as 
a decoration by atoms of the canonical cell tilings, and 
fitted to a combination of diffraction data and ab initio 
relaxations. 

Ti-Zr-Ni crystalline phase diagram - To validate 



the experimental findings as well as the diffraction- 
based structure model [js), we computed the energies of 
quasicrystal-like and competing structures and from this 
constructed the complete ground state ternary phase di- 
agram. This is the first ab initio investigation of the 
Ti-Zr-Ni alloy system. 

Our ab initio total energy calculations were performed 
with VASP @, which is a density functional code using 
a plane-wave basis and ultrasoft Vanderbilt type pseu- 
dopotentials ||. The calculations were performed us- 
ing the generalized gradient approximation by Perdew 
and Wang |)| and a plane-wave kinetic-energy cutoff of 
^cut = 302.0 eV was chosen to ensure convergence of 
the energy. The pseudopotentials for Ti and Zr de- 
scribe the 3p and 4p states, respectively, as semi-core 
states. This was found to be necessary to avoid unphys- 
ical short distances between Ti and Zr atoms. Atomic- 
level forces were calculated and relaxations with a con- 
jugate gradient method were performed. The positions 
of the atoms, as well as the shape and volume of the 
unit cells, were relaxed until the total electronic energy 
changed by less than 1 meV. This corresponds to atomic- 
level forces F max < 0.02 eV/A. The size of the k-point 
mesh was chosen to give the same accuracy for the en- 
ergy. 

The ground state phase diagram is constructed by cal- 
culating the ground state energy surface. This corre- 
sponds to determining the convex hull of the set of energy 
points as a function of the composition, as determined 
by the energy calculations for the different phases. Be- 
fore investigating the whole ternary ground state phase 
diagram, the phase diagram for the binary phases is de- 
termined from the energy calculations. For each con- 
centration the lowest energy structure is determined (see 
Table |). The convex hull of these points is then con- 
structed by considering chemical equilibria between the 
different phases. 

For the Ti-Ni phase diagram, we find in agreement with 
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TABLE I: Structures and electronic energies of the competing 
Ti-Zr-Ni phases: The calculated energies, E, the energy dif- 
ference to the groundstate, AE, and the heats of formation, 
AH{, are given. Wherever available, experimental values are 
given in parentheses p0(|. 
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the experimental phase diagram JTlj] the ground state 
phases a-Ti (hep), 7-Ni (fee), <5-Ti 2 Ni (E9 3 ), and TiNi 3 
(D0 24 ). The TiNi martensite structure B19' is found 
to be lower in energy than the high-temperature phase 
B2 and nearly degenerate with the Bf structure, which 
to our knowledge has not been observed. 

The ground state phases in the computed Zr-Ni phase 
diagram are, in addition to a-Zr (hep) and 7-Ni (fee), the 
binary phases Zr 2 Ni (C16), ZrNi (B f ), ZrNi 3 (DO19) and 
ZrNi5 (C15b). This is consistent with the binary phase 
diagram The calculated energies for both Z^Nim 

phases, Pbca and Aba2, indicate that these phases are 
slightly unstable (AE = 5meV/atom) and will decom- 
pose into ZrNi and ZrNi 3 . Our ab initio calculation of 
the binary ground state phase diagram and of the heats 
of formation (Table |jj) show excellent agreement with ex- 
periment. Furthermore, our calculated lattice parame- 
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FIG. 1: Ground state phase diagram of Ti-Zr-Ni showing 
the creases (solid lines) of the minimum energy surface of 
the competing phases. Each crease represents a domain of 
coexistence of the two phases at that crease's end points. 
Each three-phase coexistence domain is a triangle bounded 
by creases with endpoints at the corresponding single-phase 
points. Shaded regions represent the experimental composi- 
tion of the quasicrystal and the approximant phase. 



ters (not shown here) are within 2% of the experimental 
values for all observed structures. 

Besides the approximant phase W-TiZrNi and the qua- 
sicrystal z-TiZrNi, two ternary TiZrNi phases are ob- 
served for Ni concentrations below 50%, the hexago- 

[3 and the cubic <5-(Ti,Zr) 2 Ni 



nal A-TiZrNi phase (13 



phase [[l5| |l6| . The 6 phase is stabilized by small amounts 
of oxygen and is not found in samples low in oxygen [fl6[ ; 
our calculations (Table ||) found it to be unstable with 
respect to the binary phases. The atomically disordered 
C14 Laves phase A-TiZrNi is often seen in quasicrystal 
samples, but is stable only at high temperatures [0. To- 
tal energy calculations for nine deterministic ternary vari- 
ants of the C14 phase showed all to be unstable against 
decomposition into binary phases, except for the nickel- 
poor A-Ti6Zr 4 Ni 2 ; but that in turn is computed to be 
less stable than the W-TiZrNi "approximant" structure, 
which has practically the same composition E3 . 

Quasicrystal phase stability - It is nontrivial to set up a 
computation of the quasicrystal total energy, since there 
is no finite unit cell. We take advantage of the description 
of the quasicrystal as a space-filling tiling of a few types 
of cell, each type containing a fixed placement ( "decora- 
tion") by atoms. A periodic packing of the same tiles, 
after decoration, forms an "approximant" structure. To 
estimate the total energy of the quasicrystal, we calcu- 
lated that of several approximants, each made from es- 
sentially one tile type. 
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TABLE II: Electronic energies of periodic tiling structures 
approximating i-TiZrNi. Energies, E, energy differences to 
the groundstate, AE, and to the binary phases, A_Ebinary, are 
given. 
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The first, simpler tiling model |l|] uses the "Ammann" 
cells, a prolate (PR) and an oblate (OR) rhombohedron 
with edges of length 0.516 nm, as well as a composite tile 
called the rhombic dodecahedron (RD). This model's 
decoration is Ni on vertices, Ti on edges, and Zr (in the 
role of a larger atom) in interiors. Our Ammann-cell 
approximants do not contain the Bergman cluster and 
have 4 to 19 atoms per periodic unit cell. 

The second tiling model |[7| uses larger "canonical 
cells" known as A, B, C, and D. Canonical cell structures 
can be viewed as a particular way to group Ammann 
tiles, so as to maximize the frequency of the icosahedrally 
symmetric 45-atom "Bergman cluster" . The cluster cen- 
ters are the canonical cell tile corners, linked by edges of 
1.24 nm. This model's atomic decoration is locally simi- 
lar to the first model's, but permits many more site types 
depending on the local context. In a previous paper 
we determined the structure of i-TiZrNi by a constrained 
least-squares fit of this model, in which the chemical site 
occupations were refined using X-ray and neutron diffrac- 
tion data, while the atomic positions were optimized by 
ab initio relaxations of the resulting atomic decoration for 
periodic tilings; two iterations of this two-step procedure 
gave convergence. Several sites were found to have mixed 
occupation, which must be chosen one way or the other in 
the ab initio calculations, producing the numbered vari- 
ant decorations in Table ||. The canonical-cell approxi- 
mants (Aq, B2C2 and D2 packings) are relatively large, 
with 81 to 123 atoms per crystallographic unit cell || \n\ . 

From the energies of the ternary approximant struc- 
tures the groundstate energy surface is constructed in the 
same way as done for the binary phases. Table || gives 
the energy of formation AE we found for these struc- 
tures, which is defined as the energy difference to the co- 
existing mixture of three competing binary phases with 



the same composition and the lowest possible energy, i.e. 
the groundstate energy surface. Aginary is defined in 
comparison to a coexisting mixture of three competing 
binary phases, which in our calculation are a-Ti, Zi^Ni, 
and (depending on the composition) either a-Zr or ZrNi. 
Experimentally, the competing phases of the quasicrystal 
are a- Ti/Zr, Ti 2 Ni, Zr 2 Ni, A-TiZrNi (Laves phase), and 
the W-TiZrNi approximant phase Q. 

All three periodic Ammann tiling structures are un- 
stable against the competing binary phases. The OR 
packing is particularly unfavorable energetically. On the 
other hand, the PR packing has AE 0. This struc- 
ture is equivalent to the cubic C15 (MgCu2-type) Laves 
phase |lj] ; it can be obtained by replacing most of the Ni 
atoms in the C15 ZrNi 2 structure by Ti. 

The quasicrystal-like periodic canonical cell tilings, on 
the other hand, have significantly lower energies, by up to 
28meV/atom compared to the competing binary phases 
and also lower in energy than the ternary C14 phase, A- 
TigZr4Ni2, as mentioned earlier. We found the three Aq 
structures to be groundstate structures and the B2C2 as 
well as the D2 structures to lie slightly above the ground- 
state energy surface by 4 to 9 meV/atom. The A§ tiling 
corresponds to the experimentally observed VF-TiZrNi 
phase (which inspired the structure model |^]). Diffrac- 
tion experiments show that H^-TiZrNi has Ti/Zr disorder 
on the two kinds of site which lie on the plane bisecting 
the (100) axis linking two Bergman clusters there- 
fore we investigated variant chemical occupations of these 
so-called "glue" sites. We found that all three variants 
give the lowest energy for their composition. In fact, 
the three structures form a line on the groundstate en- 
ergy surface, indicating that Ti/Zr disorder on these sites 
costs no energy. 

The next larger canonical cell tilings, the B2C2 and the 
D2 tilings, do not correspond to experimentally observed 
phases, but are computed to be stable with respect to the 
binary phases. We investigated variant site occupations 
in both these structures (l9) . Changing Zr atoms to Ti on 
certain sites only weakly influences the energy, whereas 
changing Ti atoms to Ni leads to significant increase of 
the energy, indicating that the quasicrystal might contain 
Ti/Zr disorder but not Ti/Ni disorder, in agreement with 
the structural refinement of i-TiZrNi of Ref. ||. 

The stability of the decoration model of the three peri- 
odic canonical cell tilings is a strong indication that the 
decoration of a larger tiling should be practically sta- 
ble. We estimated the energy of i-TiZrNi by summing 
the energy of each constituent tile, as found from the 
structures with only one kind of tile. In an infinite icosa- 
hcdral canonical-cell tiling, the number ratio of tiles is 
N(A) : N(BC) : N(D) = (3[3 - \/5] : 1 : [y/E - 2]), 
where we adopted the "magic" value C = 3(1 - 2/\/E), 
a parameter for the frequency of D cells JlT| . Taking 
the lowest energy tiling structures Aq(1), -E^C^l) and 
the I?2(l) yields for the composition of the quasicrystal 
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Ti46.oZr37.6Nii6.4 (1), 4% higher in Ti than the exper- 
imental composition of Ti^.sZr^.sNiiy. Taking the Zr 
rich B 2 C 2 (2) structure instead yields Ti4 2 .7Zr4 .gNi 16 .4 
(2) within 1% of the experimental composition. The en- 
ergies of the periodic tilings are not enough to deter- 
mine the energy of the quasicrystal since inter-tile inter- 
actions (called "tile Hamiltonian" ) need to be considered 
too Q . Neglecting any tile-tile interactions the two qua- 
sicrystal models are 21 and 17 meV/atom lower than the 
competing binary phases. Both structures, however, are 
slightly above the ground state energy surface by 3 and 
6 meV/atom respectively pj[ . This combined with the 
experimental results suggests that the icosahedral TiZrNi 
quasicrystal really may be a ground state quasicrystal. 

What mechanism stabilizes the quasicrystal-like ap- 
proximant structures? Ti and Zr exhibit a zero heat of 
mixing. Ti and Ni as well as Zr and Ni, on the other 
hand, show a strongly attractive interaction. Thus, it is 
energetically favorable for the Ni atoms to be surrounded 
by Ti or Zr. This is reflected in a large charge transfer 
from the Ti/Zr to Ni and a strong hybridization between 
the Ni and the Ti/Zr subbands which we observe in our 
calculations. This explains why in the quasicrystalline 
structure the Ni atoms occupy the sites along the edges 
of the Ammann tiles surrounded by Ti and Zr. Further- 
more, since Zr is slightly larger than Ti it is no surprise 
that Zr occupies the more open sites of the structure. 

However, the unstable small-Ammann-tiling approxi- 
mants PR, OR, and RD have virtually the same local 
atomic structure as the canonical-cell approximants A§, 
B2C2, and Z?2- The main difference is that complete 
Bergman clusters are absent in the small Ammann ap- 
proximants, but dense in the canonical cell approximants 
(where they contain ~ 1/2 of the atoms); we conjecture 
that this is responsible for the energy difference. 

Conclusion - The first systematic investigation of the 
energies of structures of the binary Ti-Ni and Zr-Ni phase 
diagrams was presented and it was shown that ab initio 
calculations for these systems yield excellent agreement 
for the ground state phases with experiments. The phase 
transition between the W phase and the quasicrystal ob- 
served experimentally as well as the long-time annealing 
experiments demonstrate that the quasicrystal is a sta- 
ble low temperature phase. The ab initio energies of the 
nickel-poor Ti-Zr-Ni ternary crystalline phases showed 
that only the large-cell "quasicrystal approximants" con- 
taining complete Bergman clusters were stable against 
decomposition into the binary phases. This demonstrates 
that the i-TiZrNi phase, the IF-TiZrNi phase, or a sim- 
ilar large-cell approximant is the groundstate of the sys- 
tem. The calculations (Table ||) imply a total energy dif- 
ference favoring VU-TiZrNi by the order of 100 meV per 
tile. The experimental situation indicates that i-TiZrNi 
is stabilized by energy. Conceivably the quasicrystal - 
or large approximants that mix different canonical cell 
types - could be stabilized by the cell-cell interaction 



energies, which were hitherto neglected. Energy calcu- 
lations on those larger approximants would permit the 
extraction of the tile Hamiltonian, which is needed to 
decide whether the quasicrystal is a ground state, and to 
address its long-range structure. 
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